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1. Introduction 


The detection of very high energy astrophysical neutrinos in the IceCube experiment [|I|, Q] 
opens up a new possibility for investigation of particle acceleration processes in the Universe. The 
neutrino production in cosmos is possible via the pp and py interactions and the decay chains 
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e^VgV^. The observed astrophysical flux emerges from under more steep air 


show spectrum at about 50 TeV and has a cutoff at 2 PeV. The neutrino typically carries small part 
of the primary proton energy, Ey Ri O.OSEp, so the protons with energies up to E^ax ~ 10^^ eV 
are required to explain the observations (this energy is ~ 10^^ eV/nucleon in the case of neutrino 
production by nuclei). Assuming an E^'^ power-law spectrum, the measured differential flux of 
astrophysical neutrinos is E^E{E) = 2.9 x lO^^GeV^^s^^sr^' for the sum of the three evenly 
distributed neutrino flavors. The sources of observed neutrinos are not yet identified. The detected 
37 events are scattered over the sky and do not show evident correlation with any astronomical 
objects [Q]. It seems that Galactic sources might account only for a minority of events. The detected 
astrophysical neutrinos could be produced in extragalactic sources of ultra high energy protons and 
nuclei. The discussion about potential sources of very high energy neutrinos in the light of the last 
experimental results can be found in 0, ^ where other useful references are given. 

The rare extragalactic Type Iln supernova remnants are considered in the present paper as 
sources of diffuse high energy neutrinos. It is well established that supernova remnants are ef¬ 
ficient accelerators of protons, nuclei and electrons. They are the principle sources of Galactic 
cosmic rays. The diffusive shock acceleration mechanism suggested in [Q, ^ ^ |T^ can provide 
the acceleration of protons and nuclei in the most frequent Type IIP, la, Ib/c supernova events up 
to about lO^^Z eV that allows to explain the spectrum and composition of Galactic cosmic rays 
with a proton-helium knee at 3 x 10^^ eV and the maximum energy ~ 10^^ eV where iron nuclei 
dominate, see iini]- Two orders of magnitude larger E^ax ~ 10^^ eV/nucleon is needed to explain 
the IceCube neutrino observations. It can be achieved with the Type Iln supernovae that stand out 
because of extremely dense wind of their progenitor stars with a mass loss rate 10^^ — 10^' Mq 
yr^i 0- As it will be shown below, the large kinetic energy of explosion and very high gas den¬ 
sity in the acceleration region lead to the needed energy of accelerated particles and efficiency of 
neutrino production in pp interactions. 

Diffusive shock acceleration by supernova shocks propagating in dense stellar winds was al¬ 
ready considered in |0, Q. Simple analytical estimates showed that radiowaves, gamma-rays 
and neutrinos might be observable from the nearest Type Iln supernova remnants if the efficient 
diffusive shock acceleration takes place in these objects. In the present paper we develop this idea 
further and investigate whether the flux of neutrinos produced in extragalactic Type Iln supernova 
remnants can explain the IceCube data. For this purpose we perform numerical modeling of par¬ 
ticle acceleration in a supernova remnant produced by Type Iln supernova explosion and calculate 
neutrino production. Our model of nonlinear diffusive shock acceleration describes the remnant 
evolution and the production of energetic particles. The detailed description of the model was pre¬ 
sented in [ lA I and the simplified version of the model was used in [ |ll| ] for the explanation of energy 
spectrum and composition of Galactic cosmic rays. Similar numerical models of diffusive shock 
acceleration in supernova remnants were developed and employed in [ |T^ [T^ p^ . 

The paper is organized as follows. In the next Section 2 we describe modeling of particle 
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acceleration and calculate the spectrum of neutrinos produced in Type Iln supernova remnants. 
These results are used in Section 3 for calculation of diffuse neutrino background in the expanding 
Universe. The discussion of results and conclusions are given in Sections 4 and 5. 


2. Modeling of diffusive shock acceleration in the remnant of Type Iln supernova 


Details of our model of nonlinear diffusive shock acceleration can be found in []15|]. The 
model contains coupled spherically symmetric hydrodynamic equations and the transport equations 
for energetic protons, ions and electrons. The forward and reverse shocks are included in the 
consideration. 

The blast wave produced by Type Iln supernova explosion propagates through the wind of 
the presupernova star. We assume that the initial stellar wind density profile is described by the 
following expression: 


P = 


M 


47rM„ 


k2- 


(2.1) 


Here M is the mass-loss rate of the supernova progenitor star, is the wind velocity and r is the 
distance from the center of explosion. 

We use the following pai'ameters of the supernova explosion. The explosion energy Esn = 10^^ 
erg, ejecta mass Mej = 10 Mq, M = 10^^ Mq yr^', the parameter of ejecta velocity distribution 
k = 9 (this parameter describes the power-law density profile oc of the outer part of the 
ejecta that freely expands after supernova explosion), the stellar wind velocity = 100 km s^^ 
that are the characteristic values for Type Iln supernova [|T^]. 

The injection efficiency of thermal protons t] is taken to be independent of time T] = 0.01. Pro¬ 
tons of mass m are injected at the forward and reverse shocks. The high injection efficiency results 
in the significant shock modification already at early stage of the supernova remnant expansion. 

The secondary electrons and positrons are effectively produced in the dense medium of Type 
Iln supernova remnant via pp interactions. That is why we do not take into account injection of 
thermal electrons at the shocks in the present calculations. 

Cosmic ray diffusion is determined by particle scattering on magnetic inhomogeneities. The 
cosmic ray streaming instability increases the level of magnetohydrodynamic turbulence in the 
shock vicinity and even significantly amplify the absolute value of magnetic field in young 
supernova remnants [1^, It decreases the diffusion coefficient and increases the maximum 


energy of accelerated particles. The results of continuing theoretical study of this effect can be 
found in review papers |^, ^|. 

In our calculations below, we use the Bohm value of the diffusion coefficient Db = pvc/'iqB, 
where q is the electric charge of particles. The amplified upsfream magnetic field is B = Vf^/4^^p /Ma- 
Here Vf is fhe speed of fhe forward shock. The value of the parameter Ma = 10 describing the mag¬ 
netic field amplification is assumed. 

Figures (l)-(2) illustrate the results of our numerical calculations. 

We stop our calculations at ? = 30 yr when the stellar wind mass swept up by the forward 
shock is of the order of 20 Mq. This value of the total mass loss is expected if the initial mass 
of the Type Iln supernova progenitor is 60 — VOM©. The rest goes into ejecta (~ lOM©) and the 
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Figure 1: Spectra of particles produced in the supernova remnant during 30 yr after explosion. Spectrum of 
protons (thick solid line ), spectrum of secondary electrons (multiplied on 10^, thin solid line), spectrum of 
neutrinos (thick dashed line) are shown. 



Figure 2: Dependencies on time of fluxes from supernova remnant at distance 1 Mpc. We show the neutrino 
flux at 1 PeV produced via pion decay (thick dotted line), gamma-ray flux at 1 GeV produced via pion decay 
(thick dashed line). The evolution of synchrotron gamma-ray flux at 1 MeV (thick solid line) and radio-fluxes 
at 1.4 GHz, 5 GHz and 15 GHz (thin solid, dashed and dotted lines respectively) are also shown. 


black hole (~ 40Mq). At later times the forward shock enters rarefied medium created by the fast 
tenuous wind of supernova progenitor at the main sequence stage. So the production of neutrinos 
becomes negligible at this time. 

The spectra of particles produced during 30 years after explosion are shown in Fig.l. They 
are calculated as the sum of the spectra integrated throughout simulation domain and of the time- 
integrated outward diffusive flux at the simulation boundary at r = 2Rf where R / is the forward 
shock radius. The maximum energy of accelerated protons is about 30 PeV, while higher energy 
particles have already left the remnant. We also show the time-integrated spectrum of neutrinos 
produced in pp interactions. We checked that py interactions result in the negligible contribution 
to neutrino fluxes for realistic supernova optical luminosity. 

About 25% of the kinetic energy of explosion goes to cosmic rays. 
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Evolution of non-thermal emission produced in the supernova remnant at distance 1 Mpc is 
shown in Fig.2. We take into account the synchr'otron self-absorption and free-free thermal absorp¬ 
tion that are important for radio-supernovae [|^]. We use the temperature T = lO^K of circumstel- 
lar wind that is a characteristic value for dense stellar winds ionized by radiation coming from the 
forward and reverse shocks [M|. 

It is instructive to compare results of the present calculations with the approximate analytical 
expression for the maximum energy of accelerated protons derived in our paper [|n|]. The maxi¬ 
mum energy of protons Emax accelerated by a supernova shock in the stellar wind was estimated 
as 


= 8PeV 


M 


10-4 Mq yr- 


1/2 


Mvi 


100 km s“ 


-1/2 


f E: 


/ ^SN 
\ 10^2 erg 


( Mej 
VIOMq 


( 2 . 2 ) 


A similar estimate for Type la supernova explosion in the uniform medium gives the "knee" 
energy ~ 3 PeV for protons. The shocks propagating in stellar winds accelerate particles to higher 
energies. In particular Type Iln supernova remnants with their dense winds with mass-loss rate M ~ 
10^^ Mq yr-' can accelerate protons up to 80 PeV according to Eq. (2.2). This order of magnitude 
estimate is in agreement with our numerical results illustrated in Fig. 1. Thus, explanation of 
proton acceleration up to ~ 10^^ eV in Type Iln supernova remnants is in line with the acceleration 
of cosmic rays up to the knee in Type la supernova remnants. The "knee" energy of accelerated 
particles is reached at the beginning of the Sedov stage when a significant fraction of explosion 
energy transfers to the shock wave. In denser wind, the amplified magnefic field is sfronger fhaf 
leads fo larger Emax- 


3. Calculation of background neutrino flux 

The neufrino specfrum produced by a single supernova remnanf Q{Ey) (see Fig.l) can be used 
for fhe deferminafion of exfragalacfic neufrino background. Disfribufed in fhe Universe Type Iln 
supernova remnanfs give fhe following diffuse flux of neufrinos: 


E(Ev) = 


AkHq Jo 


■Z,„ax Q{[\+z)Ey 

az 


IVm(l -t-z)'” 


Va™(l+z)3 + aA 47r//o Je, 


rf+z,„ 

Jev 


v)Ev , 

dE' 


VsnQ{E') 


m+i 


\//Ey + 

(3.1) 

Here fhe adiabatic energy loss of neufrinos produced af fhe redshiffs 0 < z < Zmax is taken 
into account. The present neutrino production rate per unit energy and volume is VsnQ{Ev), where 
Vsn is the rate of Type Iln supernovae at z = 0 while the cosmological evolution of the sources 
in the comoving volume is described as (1 +z)’" (m = 0 implies no evolution). The evolution 
parameter m = 3.3 for z < 1 and no evolution at z > 1, the maximum redshift Zmax = 5 and the 
rate V™ = 10^^ Mpc^^ yr-^ at z = 0 are assumed in our calculations (see e.g. 125]). This rate of 
Type Iln supernovae is 100 times lower than the rate of all core collapse supernovae. Hq = 70 km 
s^^ Mpc^' is the Hubble parameter at the present epoch, the matter density in the flat Universe is 
Q-m = 0.28, and the A -term is Q.\ = 0.72. 
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Figure 3: Calculated spectra of neutrino produced by Iln SNRs in the expanding Universe (solid line). 
IceCube data Q are also shown. 


The calculated background neutrino spectrum is shown in Fig.3. The figure demonstrates very 
good fit of our calculations to the IceCube data. 

Using the same approach we can calculate the flux of exfragalacfic profons. Our resulfs are 
compared wifh cosmic ray dafa in Fig.4. So fhe profon flux produced by exfragalacfic Iln supernova 
is below fhe measured cosmic ray fluxes. 

The possible confribufion of Galacfic Type Iln SN evenfs fo fhe observed cosmic ray infensify 
mighf be comparable wifh fhe Type Ilb supernovae as if was calculafed in our paper [11]. Bofh 
types of supernovae probably produce approximafely fhe same amounf of cosmic rays wifh almost 
the same maximum energy of particles but SNIIn generate much more secondary particles and 
neutrinos. The low expected rate of these Galactic supernovae, about one per 5 thousand years, 
makes estimates of their contribution to the observed intensity of very high energy Galactic cosmic 
rays rather uncertain. 

Simple order of magnitude estimates can be done to clarify how the obtained neutrino flux 
depends on supernova paramefers. The main production of high energy particles and neufrinos is 
finished af fhe beginning of the Sedov stage when the shock radius Rg can be determined from 
the condition Mgj = 4k drr^p = MRs/uw The spectrum is assumed for the accelerated 
particles. The neutrino energy flux expected from a supernova at distance D can be estimated as 
(see also [ 0 , 111 ) 


2 cVfOppM 

= 16712 ln(F^ax/'Hc2) mulRfD^ 


Mej 

lOMoy 

(3.2) 

Here ^CR is the ratio of cosmic ray pressure to the ram pressure of the shock pVp Ky « 0.25 is 
the fraction of energy that goes into neutrinos in pp interactions, Opp is the cross-section of pp 
interactions and Fmax is the maximum energy of accelerated protons given by Eq. (2.2). The 
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expression is valid when the shock radius Rf > /?min = cGppM/4nu„rnVf and the energy losses via 
pp interactions are smaller than the adiabatic losses. 

Expression (3.2) can be used for the determination of neutrino background produced by all 
Type Iln supernovae: 


P(j7\j72_ ^^CrKv VsnC^VfOppM^ . Rs 

'' 16712 ln(£'n,a^/mc2) Homul ^ Rmin 


10 ^'ergcm ^s ^sr 'x 


£ f ^ \ ( uw \ Y esn 

^0 yr^ V Yo^o/ YoOkms-V Yo^^ erg^ Yo-^Mpc-V^^ 

(3.3) 

The background flux is higher if the cosmological evolution of supernova rate is taken into account 
(see Eq. (3.1)). These estimates are valid at neutrino energies Ey < O.OSEmax- They confirm 
the efficiency of high energy neutrino production in supernova explosion in very dense wind of a 
progenitor star, as it takes place in Type Iln supernovae. 



4. Discussion 


Our calculations show that the IceCube data can be explained by the neutrinos from Type 
Iln supernovae. If so, the arrival direction of every IceCube neutrino should coincide with the 
direction to some Type Iln supernova. However, we observe Type Iln supernovae with the redshifts 
z < 0.1 while the background is determined by the supernovae with z ~ 1. Therefore we expect 
that only several percent of observed neutrinos i.e. about one IceCube neutrino event are associated 
with the known Type Iln supernovae. Unfortunately, the anival directions of majority of IceCube 
neutrinos are known with low angular resolution (about 15 degrees). We found that several Type Iln 
supernovae were indeed observed in the vicinity of some IceCube neutrinos. However this result 
has low statistical significance. 

Nevertheless we note that the arrival direction of one IceCube PeV neutrino candidate (event 
20) is within 5 angular degrees from the SN 1978K. This nearest Iln supernova at distance 4 Mpc 
is for decades observed in radio. X-rays and optics (e.g. [^]). The mass-loss rate of SN 1978K 
progenitor is estimated as M = 2 • 10^^ Mq yr^' that is 5 times lower than we use in our 
calculations. Thus the fluxes expected from this supernova are not so high as ones shown in Eig.2. 
In spite of this it is possible that the PeV neutrino candidate was emitted by SN 1978K. 

The arrival direction of another IceCube PeV neutrino candidate (event 35) is within 10 angular 
degrees from the SN 1996cr. This nearest Iln supernova at distance 4 Mpc was also detected in 


radio. X-rays and optics (e.g. [|29|]). It is known that the forward shock of SN1996cr entered the 
dense shell of circumstellar matter. The radio flux of this supernova is only several times lower 
than the flux shown in Pig.2. That is why this IceCube PeV neutrino candidate could be emitted by 
SN 1996cr. 

On the other hand, 8 of 37 IceCube neutrinos left tracks. The arrival directions of these neu¬ 
trinos are determined with good angular resolution (about 1 degree). We found that none of these 
neutrinos is associated with known Type Iln supernovae. The expected number of coincidences is 
of the order of 0.3. Thus, only one track neutrino can be observed from the vicinity of Type Iln 
supernova during one decade of IceCube operation. 
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Figure 4: Comparison of calculated spectrum of cosmic ray protons produced by extragalactic Iln SNRs in 
the expanding Universe (gray solid line) and available cosmic ray data 


5. Conclusions 

Our main conclusions are the following: 

1) The diffusive shock acceleration of particles in Type Iln supernova remnants and the pro¬ 
duction of neutrinos via pp interactions in the dense presupernova winds can explain the diffuse 
flux of high energy astrophysical neutrinos observed in the IceCube experiment. 

2) The calculated maximum energy of protons accelerated in remnants of Iln supernova is 
close to 10^^ eV. This value is higher than the maximum energy achieved in the main part of SNRs 
and is explained by the high density of the circumstellar matter. 

3) The efficient acceleration of particles and production of secondary electrons and positrons 
result in the fluxes of radiowaves and gamma-rays that can be observed from the nearest Type Iln 
supernova remnants. 

4) Future IceCube operation and the search of correlations between neutrino arrival directions 
and the directions to Type Iln supernovae may check whether these objects are the sources of high 
energy neutrinos but only one track event during 10 years is expected from these supernovae with 
the redshifts less than 0.1. The bright phase of a Type Iln supernova remnant as a source of PeV 
neutrinos lasts for about several years after the supernova explosion. 

The work was supported by the Russian Foundation for Basic Research grant 13-02-00056 
and by the Russian Federation Ministry of Science and Education contract 14.518.11.7046. 
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